* Limits
gisvol.up(t)  = 99.9;
gisvol.lo(t)  =  .1;


* Initial conditions
*GISVOL.fx(tfirst) = volzero;
